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Abstract 

Pauli-projected random gaussians are used as a representation to 
solve the shell model equations. The elements of the representation are 
chosen by a variational procedure. This scheme is particularly suited 
to describe cluster formation and cluster decay in nuclei. It overcomes 
the basis-size problem of the ordinary shell model and the technical 
difficulties of the cluster-configuration shell model. The model repro- 
duces the a-decay width of 212 Po satisfactorily. 



1 Introduction 

By complementing the shell model with a cluster-model basis, the a-decay of 
212 Po can be described very well Jj] . This cluster-configuration shell model or 
hybrid model is, however, difficult to apply. We have therefore developed a 
simpler model which is closer to the conventional shell model but in practical 
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applications may prove to be far superior. The objective of this paper is to 
report on this new model and its performance in the calculation of the ground 
state of 212 Po. To this end, we compare its results with those of the hybrid 
model 0. 

The new model is essentially a shell model, in which the valence orbits are 
replaced by combinations of gaussian functions of different size parameters, 
with exact Pauli projection off the core orbits. These gaussians are more 
flexible than the harmonic oscillator single-particle (s.p.) functions because 
their width parameters are free. A variational backgroud can guarantee that 
they will be automatically chosen so as to suit to the problem considered. In 
particular, the gaussians can produce a more realistic fall-off in the surface, 
which is important for the description of decay processes. Moreover, in the 
interior the gaussians are able to simulate the effects of several oscillator 
shells. Note that the Pauli projection provides the gaussians with some, 
say n , nodes, so that n gaussians in the expansion of a s.p. function will 
involve s.p. orbits of no, . . . , no + n — 1 nodes. In spite of the inclusion of 
such high-lying shells, the problem may still remain tractable if the basis 
elements are carefully selected. We admit states in the basis depending on 
their contribution to the g.s. energy. This method of choosing random states 
as candidates for the basis is called || "stochastic variational method" . 



2 The formalism 

The basic building block of our s.p. basis is the gaussian (or nodeless har- 
monic oscillator function) of size parameter z/, defined by 

As any well-behaved function can be approximated with any prescribed 
precision by a linear combination of gaussians of different size parameters, 
X)i c iVT{ r )i one can use sucn combinations to approximate the valence s.p. 
functions of the shell-model Hamiltonian. 

These gaussians are not orthogonal to the core orbits. To take care of 
antisymmetrization, one can introduce the operator that projects out the s.p. 
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states occupied in the core [|: 



j=5 



(2) 



The index i distinguishes the proton {i = 1,2) and neutron {i = 3,4) orbits. 
The Pauli-correct s.p. wave functions are then constructed as 



j rn 



(3) 



where is the spin function of the zth nucleon. 

In the description of 212 Po it is reasonable to assume the core, 208 Pb, to 
be passive and massive. (This approximation could be removed by using the 
formalism of the cluster-orbital shell model M.) 

As these s.p. functions are fully analogous to those in the conventional 
shell model, the many-particle basis can be constructed in the usual way. 
In proton-neutron formalism and jj coupling scheme the two-proton and 
two-neutron functions are given by 
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A 



and 



(4) 



(5) 



where the operator A antisymmetrizes between like particles, 7r = {h,ji, h, J2, J-*}, 
v — {h,j3, U,j4, Ju}- A four-particle state is formed as 



{■kv)JM 



7T V 



JM 



The trial function is then 



JM 



(ttu)JM- 



(6) 



(7) 
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To describe the dynamics of the valence particles, we minimize the ex- 
pectation value of their Hamiltonian 



h = j2m + u t )+ £ % 

i=i i<j<i<4 



(8) 
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where T is the kinetic energy, U is the core-particle interaction and V is the 
particle-particle interaction. 

The trial function ^ jm contains a great number of terms belonging to 
different configurations and size parameters. Since the dependence of the 
eigenvalue problem upon the variational parameters q is linear one can de- 
termine them by finding the lowest eigenvalue of a generalized linear algebraic 
eigenvalue problem, i. e. 

He = E N Ac, (9) 

where H and A are the Hamiltonian and overlap matrices. The subscript iV 
stands for the dimension of the eigenvalue problem. The overlap matrix A 
is nondiagonal due to the nonorthogonality of the functions f[„) JM . 

This variational approach is only expected to yield a good approximation 
to the g.s. energy if the size parameters v\,...,v\ are adequately chosen. 
Previous experience with the stochastic multicluster model shows that 
the number of basis states needed is not excessive even in the five-cluster 
case, which is comparable with the present core plus four-nucleon system. 

We generated trial basis states randomly, and singled out, by an admit- 
tance test, the ones whose contribution to the g.s. energy was larger than 
a preset limit, e ("utility testing" 0). That is, we start (k=l) by making 
a random choice of the size parameter quartet v\ , . . . , z/| from the physically 
important v region, and set up a basis state of a particular configuration. 
The energy expectation value gives the first approximate energy, E\. Then 
we generate a new random set uf, vf, and solve eq. (§) for N — 2. Due to 
the variational character of the method any parameter quartet improves the 
energy (E 2 < Ei), but we only adopt the basis element if E\ — E 2 > e. This 
procedure is then repeated with new random elements up to convergence, 
and then the basis is enlarged by including further configurations. 

At each trial step the energy can be determined by solving the eigenvalue 
problem using any conventional numerical algorithm. But as we are only 
interested in the g.s. energy, and we proceed stepwise, we can use the much 
more economical modified Jacobi method PJ. The Jacobi method (see e. g. 
Ref. 0) solves the eigenvalue problem HC = EC eliminating the nondiag- 
onal elements by successive plane rotations. In the modified Jacobi method, 
the lowest eigenvalue is only determined by elimination of the nondiagonal 
elements corresponding to that eigenvalue ||. Thus it uses, in every step, 
the results of the diagonalization of the previous step. Although the modi- 
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fied Jacobi method can be generalized to solve fully the eigenvalue problem 
|J, we found it more convenient to re-orthogonalize our basis in each step 
of the random search procedure. This is advantageous because A is block- 
diagonal (each configuration provides a block) and we need to orthogonalize 
only within the configuration whose basis is being enlarged. Since the number 
of basis states in any configuration is low, this re-orthogonalization requires 
almost no computational effort. 

A delicate point of this procedure is the choice of e because a too large 
value bogs down the procedure, whereas a too small value results in highly 
redundant bases of large dimensions. It is expedient to re-adapt the 
acceptance criterion dynamically. For instance, in the present calculation we 
started with e = 0.01 MeV, and after 10 failed attempts to find a good enough 
basis state, we reduced e by 2, and this was repeated two more times. In this 
way we managed to find further states of some significance. This procedure 
has been thorougly tested [Q] , and it was found that it is extremely powerful 
and reliable f2j. 

3 Details 

The s.p. potentials include scalar and spin-orbit Woods-Saxon parts as well 
as the Coulomb interaction. To facilitate the analytical calculation of the 
matrix elements, all potential terms are expanded in terms of gaussians. To 
be able to compare the "pure shell model" [TLJ with the present model, we use 
the potential parameter set B of ref. |1[] and employ the same Pauli projec- 
tor as there, that is the occupied s.p. states are represented by oscillator 
functions of size parameter vq = 0.083fm -2 . 

For the valence orbits we selected the size parameters from the physically 
important interval 4 fm < v~ x l 2 < 14 fm. The following s.p. quantum 
numbers were used: 

1,3 = hg/2,f7/2,iia/2,f 5/2, Pap (protons) 

= 99/2, hi/2, hh/2, 4/2, si/2,£77/2 (neutrons). 

Just as in the former model, we set up the four-particle basis by coupling all 
possible two-particle combinations up to J n = 8 and J v = 8. 

We build up the basis by the method skecthed above in each config- 
uration {7rz/}. We ordered the configurations according to their expected 
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contributions to the energy, i.e. according to the sums of the s.p. energies 
involved and, among states involving the same s.p. states, to the two-particle 
angular momenta in increasing order. Thus we begin with the configura- 
tions [[h 9/2 } 2 j7i [g 9/2 } 2 jJo, with J w = J v = 0,1,2,..., [[/tyaljjiu/aljjo, with 
■h = Ju = 0, 1, 2, . . ., [[/7/a]j w [fffl/2]jJo, with J v = J v = 0,1,2, .. . etc. If a 
configuration cannot contribute to the g.s. energy by at least e/2 3 , then it is 
omitted. 



4 Results 

First we solved the eigenvalue problem of the s.p. part of the Hamiltonian 
on the basis (|3|), The scalar depth was adjusted so as to best reproduce 
the experimental s.p. energies, as done in ref. We obtained this with 
the scalar strengths V p = 61.60 MeV and V n = 43.80 MeV for protons and 
neutrons, respectively. The resulting s.p. energies are in excellent agreement 
with the ones obtained by a direct numerical integration. 

We then calculated the g.s. (J = 0) energy of 212 Po. We terminated the 
random search at a basis dimension N = 397 and a number of configurations 
n — 118. As the number of basis states is increased the g.s. energy converges 
to E — —19.25 MeV (the experimental value is —19.35 MeV.), as seen in 
figure 1. This may be compared with the value E = —18.96 MeV of the pure 
shell model or the value E = —19.18 MeV in the hybrid model 0. 

To check the convergence we repeated the calculation using different ran- 
dom parameter sets with dimensions N = 413, n = 120 and N = 380, 
n = 127. The three energies agree within 0.04 MeV. To see the role of the 
order of the configurations, we repeated again the calculations with reversed 
order. The resulting energies are very close to the previous ones, i. e. —19.24 
MeV, —19.28 MeV and —19.25 MeV, but the basis sizes are almost doubled. 
To test the choice of e, we made three calculations with e = 0.005 MeV. 
The lowest calculated energy obtained is E = —19.27 MeV, in good agree- 
ment with the former calculations, but the basis sizes increased considerably 
( N = 631, n = 193). 

We also calculated the spectroscopic (or formation) amplitude defined by 

0(r) = ($ w |*oo), (10) 
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with 

$ (r) = $ / (r ~ rca V LM (f co ), (ii) 

' ca 

where $ Q is the intrinsic wave function of the a particle constructed from 
a Os harmonic-oscillator Slater determinant with size parameter v a . (See 
Ref. fl|] for details.) The amplitudes of the present model and of the cluster- 
configuration shell model are compared in Fig. 2. The main peak is smaller 
and its position is somewhat closer to the origin (at 8.0 fm) in the present 
model than in the cluster-configuration model (8.2 fm). 

The present spectroscopic factor is 0.011, while it is 0.023 in the hybrid 
model. 

To calculate the decay width T, we use the i?-matrix formula || 

r = 2P ^2K~a^ (12 > 

where M Q is the reduced mass, a is the channel radius and P is the Coulomb 
penetration factor. The width, calculated in the interval lOfm < a < 12 
fm, where the effect of the nuclear forces and Pauli exchanges are negligible, 
is T = 7.2 x 10 -16 MeV. This value is about one third of the hybrid-model 
value (r = 2.1 x 10~ 15 ) and about one half of the experimental (r = 1.5 x 
10~ 15 ) value. Note that the conventional shell model, like that of ref. 
give much smaller and spacially confined amplitudes, while the shell models 
with sophisticated truncation schemes |I0, [UJ give amplitudes of comparable 
extensions but smaller spectroscopic factors and decay widths. 



5 Summary 

With the help of a stochastic variational method, a flexible shell-model wave 
function has been constructed for the g.s. of 212 Po using combinations of 
gaussian functions as s.p. orbits. The present model gives nearly as good a 
result for the a-decay width as our former cluster-configuration shell model, 
which is sharpened just to describe a-decay. The chief merit of the model 
presented here is its technical simplicity, which makes it easily applicable to 
other a-decays or even to heavy-cluster decays. 

This work was supported by the OTKA grants (No. 3010 and F4348). 
We are grateful to Prof. R. G. Lovas for useful discussions. 
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Figure captions 



Figure 1 Convergence of the energy of 212 Po (g.s.) as a function of the basis 
dimension N. 

Figure 2 Formation amplitude g(r) (eq. (10)) as a function of the distance 
between the decaying a-particle and the daughter nucleus 208 Pb (g.s.) 
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